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Abstract 

The life-times due to Auger-electron emission for a hole on a deep elec- 
tronic shell of neutral and charged sodium clusters are studied for different 
sizes. We consider spherical clusters and calculate the Auger-transition prob- 
abilities using the energy levels and wave functions calculated in the Local- 
Density- Approximation (LDA). We obtain that Auger emission processes are 
energetically not allowed for neutral and positively charged sodium clusters. 
In general, the Auger probabilities in small Na^ clusters are remarkably dif- 
ferent from the atomic ones and exhibit a rich size dependence. The Auger 
decay times of most of the cluster sizes studied are orders of magnitude larger 

than in atoms and might be comparable with typical fragmentation times. 
36.40.-c, 32.80.Hd 
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I. INTRODUCTION 



The decay of an electronically excited cluster may take place by emission of electrons, 
atoms or photons. It is generally accepted that the radiative cooling is the decay channel with 
the largest time scale p]. Another deexcitation channel is the evaporation of atoms, which 
occurs within a characteristic time which ranges between pico- and milliseconds, depending 
on the bonding character of the clusters and the excitation energy. The fastest relaxation 
channel is the emission of photoelectrons. In addition, and like in atoms, secondary electrons 
can be emitted if the excitation energy is larger than the ionization potential of the cluster. 
These electrons result from intraband Auger processes, i.e., Auger decay of a valence hole 
and subsequent emission of a valence electron. 

In excited light atoms, the nonradiative decay through emission of Auger electrons has a 
much larger probability than the radiative decay and dominates the relaxation process 0. In 
clusters there are so far neither experimental nor theoretical studies of the intraband Auger 
probabilities. Recently, two important experimental studies on the deexcitation channels of 
small clusters after optical excitation have been performed, from which the contribution of 
Auger-processes can be inferred. Gantefor et al |J analyzed the kinetic energy distribution of 
electrons emitted from optically excited clusters. They found that the photoelectron spectra 
show contributions from three different processes: direct emission, thermionic emission, 
and "inelastic scattered electrons". The latter appear within a definite and narrow range 
of kinetic energies. One could in principle interpret the inelastic scattered electrons as 
coming from intraband Auger-processes. Reiners and Haberland || studied the competition 
between electron and atom emission after photoabsorption as a function of the photon energy. 
They observed that atom emission occurs within an energy range which is smaller than the 
cluster bandwidth. This means that excitations consisting of a hole in the bottom of the 
valence band (i.e., in the deepest electronic shells) do not lead to atom emission, despite 
that their energy is larger than the binding energy. The authors conjectured that this is due 
to intraband Auger-processes which lead to a very short life-time of such holes. 
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In this paper we present the first calculations of the intraband Auger time scales in 
metal clusters and show that Auger emission probabilities strongly depend on size and are, 
in general, orders of magnitude smaller than atomic Auger probabilities. 



II. THEORY 

We consider only closed- shell Na^v, Na^ and Na^ clusters and make use of their spherical 
symmetry. For the description of the electronic structure of the clusters we use the jellium 
model ||. Through the Auger process, a vacancy in a state \n"£"), with £ n »^» < ep (s,f 
being the Fermi level), is filled by an electron coming from a higher bound level n£. The 
energy released by this transition is transferred to a second electron, initially in a bound 
state n'£', which is ejected (i.e., excited into a continuum state k). Thus, the initial state 
consists of one hole in a bound state, and the final state consists of two holes in bound states 
and one electron in a continuum state. One can, however, consider the Auger process as 
a two hole — ► two hole transition, where the initial state is given by = c„»^» CT Ck \^k), 
with \ty k ) = c^j^o) an d £k > 0, and the final state by \^Sff) = c n £ a c n 'i/ a i \^k), with 
Sn'i'i £ ne < £f- k denotes the continuum state of the Auger electron and n"£",n£ stand for 
the states corresponding to the two final holes, j^o) refers to the electronic ground state of 
the cluster. The Auger transition probability can be calculated using Fermi's golden rule, 
and is given by 
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where p(Ef) corresponds to the density of final states and V is the operator describing the 
Coulomb interactions, which is written as 

V = \ ^234 4 ai 4 a2 C 3a2 C 4(Tl , (2) 

" 1234 

CT 1CT2 

where the sum is over the cluster energy levels. In Eq. (fj) the quantities V1234 are the 
Coulomb matrix elements 
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^1234 = // ^*(r )V 2 *(r') Mr')Mr) dfdf', (3) 

expressed in the basis of eigenfunctions corresponding to the bound and continuum levels 
of the cluster. 

In order to calculate the Auger-emission probability of Eq. (|l|) we have first to determine 
the wave functions of the three bound states and that of the continuum state which is 
involved in. This can be done by performing an approximation which is widely used in 
atomic physics || , and consists in taking the wave functions and energy levels of the original 
system, i.e., the cluster before the creation of the initial vacancy. Therefore we solve the 
Kohn-Sham equations for the spherical metal cluster, given by 

R n e(r) = e nl R nl {r). (4) 

Here, V e g (r ) is the effective one-electron potential 

V eff (f) = vj(f) + J df> + ±E xc [n<?)] , (5) 

where Vj(f) is the electrostatic potential created by the jellium distribution of charge and 
Exc[ n {f)} the exchange and correlation term ||. n(f) is the electronic density. From 
Eqs. (fD one obtains the bound states of the cluster. The numerical algorithm to solve 
Eqs. (|j) imposes the condition of regularity at the origin R n i oc r e and fixes the number of 
nodes n of R n £. There is only one solution for a fixed n and I which vanishes exponentially 
at infinity and the algorithm iterates to find it, as well as its associated eigenvalue e n g. The 
number of occupied shells determines the number of possible Auger processes for a given 
hole in the shell n"£". Since the total energy is conserved during the transition (Ei = Ef), 
the kinetic energy of the emitted Auger electron is given by = |e n "^»| — \e n 'i'\ — \£ n t\- 
For a Auger emission to take place it must obviously hold that Eki k > 0. Since also the 
total angular momentum is conserved, the two-hole final state must have the same angular 
momentum as the two hole initial state (Lj = Lf = L). This requires for the angular 
momentum of the Auger electron the condition \L — £"\ < £k < L + £", where L must 
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satisfy the inequality \£ — £'\ < L < £ + £'. Similar constraints are fulfilled by the spin 
of the Auger electron. For the calculation of the Auger continuum wave function we also 
use Eqs. (|]). However, in this case the number of nodes is unknown and the continuum 
energy eu h and the multipolarity £k are fixed. The same regularity condition is used to start 
the integration outwards, from the origin up to a fixed large radius Rq. We normalize the 
outgoing Auger wave function ~ exp(ikr)/r within a sphere of radius Rq. Thus, the density 
of final states is given by p{Ef) = Ro/(27rhv), where v refers to the velocity of the Auger 
electron ||. 

In order to determine the Auger probability we first separate the matrix elements @ 
into radial and angular factors. This is achieved by performing the multipole expansion 



l/ri 2 = J2fj,u r v K I rlCl^MC^iMi with C, v = being Y^B,® 

spherical harmonics. Evaluation of the angular factors depends on the choice of the angular- 
momentum coupling scheme. Since for the valence electrons of clusters the spin-orbit cou- 
pling is negligible, the initial and final two-hole states of the cluster can be expressed in the 
(LSJM) representation. The total transition probability into all possible states of L and 5* 
is then given by H] 
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where the functions D u and J u are the direct and exchange radial matrix elements, and the 
angular factors d u and e u are given by 

d v = (-1)W+ L (£"\\C»\\£) (£ k \\C»\\£'} { 4 L ) (7) 

ye £ v 

and 

{£" £ L 
" \- (8) 
£ £' v 

Here, (£||C w ||i") is the reduced matrix element of the spherical harmonic, multiplied by 
[47r/(2z. + l)]V2 . 



III. RESULTS AND DISCUSSION 



We have calculated the Auger transition probabilities for spherical (closed-shell) Najv, 
Na^ and Na^ clusters with 19 < N < 253. For each cluster we have determined first the 
effective potential V e g (r ), the bound Kohn-Sham states ({e n e} < 0) and the corresponding 
radial eigenfunctions {i? n ^(r)}. For a given vacancy n"£" we calculated the number of 
possible final states {n£, n'£'}. For each of the final states we determined the energy of 
the emitted electron and the different possible values of its angular momentum £ k . Then, for 
each £ k , we determined the corresponding outgoing radial wave function. This allowed us to 
calculate the Coulomb matrix elements and, using Eq. (P), the probability w(n£,n'£ f ,n"£"). 
Finally, we have calculated the total probability for an initial vacancy n"£" as W(n"£") = 
Y^u' w{n£,n'£' ,n"£"). The total probabilities W{n"£") can, of course, also be expressed as 
an energy width [E A (n"£" ) = h/W(n"£")] or as a life-time [r A (n"£") = W{n"£")- 1 }. 

The first interesting result of our study is that, within the approximation used in Eq. (5) 
for the exchange and correlation term (LDA), for neutral and positively charged spherical 
sodium clusters intraband Auger processes are energetically forbidden. There is no possible 
transition n"£", n£, n'£' yielding e k > 0. This means that, according to our LDA calculations, 
excited spherical Na^y and Na^ clusters can only decay via fragmentation (evaporation) or 
photon emission. 

Our results for negatively charged sodium clusters indicate that, in contrast to what 
occurs for Na^- and Na^, nonradiative electronic decay is possible. Due to the presence 
of the extra electron and the consequent extra Coulomb repulsion, the binding energy of 
the electrons in negatively charged clusters is smaller than in neutral or positively charged 
ones. As a consequence the whole band of bound states is shifted upwards and makes 
possible Auger transitions with e k > 0. An example of such Auger transitions for negatively 
charged sodium clusters is illustrated schematically in Fig. 1. Note that the Kohn-Sham 
effective potential V e ff(r) shows a barrier for negatively charged clusters. This might have 
consequences for the magnitude of the nonradiative emission probabilities. For instance, the 
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wave function of an emitted Auger-electron with positive energy but smaller than the energy 
barrier could have a large weight inside the cluster due to trapping effects and influence the 
transition matrix elements used to calculate W. For holes in the first two shells of Na^r 
there are many possible transitions. For instance, there are 6 ways of filling a Is- vacancy in 
Najg by emitting an electron. For Naj^ the number of such transitions is 11; for Naj~ 37 , 14, 
and for Na^ 53 , 19. 

Our calculated Auger widths of Na^ clusters are, for some cluster sizes, of the order 
of 10 _1 eV, i.e., as large as for light atoms. In Fig. 2 the Auger life-time of a hole in 
the first (Is) and second (lp) shells of spherical Na^ clusters is shown as a function of the 
cluster size. The size dependence of 7Vi(ls) and r^(lp) is very rich and shows no monotonical 
behavior. Furthermore, the life-times oscillate over many orders of magnitude. For instance, 
an initial ls-hole in excited Na^g and Nag 7 lives only few femtoseconds, almost as short as 
a Is- vacancy in light atoms with Z < 10 @. For these clusters the Auger emission is 
faster than any other deexcitation mechanism. This means that one can separate the time 
scales for the electronic and atomic relaxation. Thus, if one is interested in studying, for 
instance, the atomic motion after optical excitation of Na^g and Nag 7 , one can assume that 
the Auger process occurs immediately after the excitation and has no further influence on 
the fragmentation (evaporation) behavior. On the other hand, r^(ls) for Na^ 7 and ta[\p) 
for Na^g and Na^~ 7 are remarkable large, of the order of nanoseconds, i.e., larger that the 
life-times for any other deexcitation channel. For this other extreme case one can again 
separate the time scales and assume that the hole lives infinitely long (compared with the 
atomic relaxation). 

However, it is not always possible to perform this separation of time scales. Fig. 2 also 
shows that for most cluster sizes ta lies between pico- and nanoseconds. This magnitudes 
are comparable to the Auger life-times of positively charged ions approaching a metal surface 
at a distance of at least 2 A 0. This time scale is probably in the range in which excited 
clusters fragment ||. Thus, results of Fig. 2 suggest that there might be competition between 
Auger-emission and fragmentation channels for a vacancy in a deep shell of a Na^ cluster. 
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In Fig. 3 the Auger life-time for initial vacancies in the different electronic shells of Naj~ 37 
are shown. The solid line shows results obtained from the calculations as described before, 
whereas the dashed curve shows results obtained including relaxation effects in the final 
state (shake-off). The final state of the cluster corresponds actually to a system with N — 1 
electrons. This should have an influence on the Kohn-Sham levels, reflecting the fact that the 
many electron system relaxes due to the excess nuclear positive charge. For the cluster sizes 
considered, which have closed shells in the initial state, the final state with N — 1 electrons 
is no longer a closed-shell system and cannot be calculated using the spherical jellium model 
||. Thus, in order to take into account the relaxation effects we performed the following 
approach, commonly used in atomic physics [fj. We considered for both the initial and 
the final state a cluster with N electrons, i.e., a closed-shell configuration. However, for 
the final state we solved the Kohn-Sham problem for a positive jellium background with 
charge Q = \e\(N + 1). In this way we simulated the excess charge. The relaxation effects 
calculated within this approach do not change the qualitative trend of the results, neither for 
the size dependence nor for the shell dependence (for fixed size) of the Auger life-times, as 
seen in Fig. 3. Regarding the shell-dependence of t a for the different cluster sizes studied, 
there is no clear dependence on the level of the initial vacancy. ta(1s) is in most cases 
the smallest life-time. For some clusters r^m) shows an alternation for increasing shell 
number m = 1,2,.., like for Na5~ 37 . For other cluster sizes there is a monotonic increase 
of r(m) with m. In general one would expect a monotonic increase of the life-time, since 
the Coulomb matrix elements [Eq. @] should decrease with increasing kinetic energy of the 
Auger electron. However, due to the potential barrier shown in Fig. 1 and discussed below, 
for some values of Ek (resonances) the corresponding wave function could have a particular 
large weight inside the cluster, giving rise to deviations from the monotonic dependence of 
t a with energy. The potential barrier for Naj~ 37 is 0.7 eV high. 

In Fig. 4 we show the distribution P(e) of the emitted Auger-electrons as a function of 
their kinetic energies for Naj~ 37 and Na^. P(z) is calculated as 
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47 w(n£,n'£',n"£") 
{£) 27r^(e-e„J a + (7/2) 2 ' U 

where eki k (n£,n'£' ,n"£") is the energy of the electron emitted in the (n£,n'£',n"£")- Auger 
process, and the width 7 is taking to be 0.05 eV. A is a normalization constant. In both 
cluster sizes P(s) is dominated by electrons originated in few transitions with large probabil- 
ity. For Na^ 37 the large peak at 1.23 eV corresponds to a (Is, 3p, 3p)-transition, whereas the 
smaller peak at 0.66 eV results from a (Id, 3p, 3p)-process. For Na^ 53 the Auger spectrum is 
dominated by the (lp, 2h, 2/z)-process with a kinetic energy of 0.5 eV, while a smaller peak 
appears at 0.99 eV which corresponds to the (Is, 4s, 2/i)-transition. These two examples are 
consistent with the intuitive idea that the most probable Auger processes are those involving 
two electrons at the Fermi-level (3p for Naj~ 37 , and 2h in the case of Na^)- Results of Fig. 4 
remain unchanged if we take into account shake-off effects. Note that the kinetic energies 
of the emitted electrons are, for all clusters studied, not smaller than 0.4 eV and not larger 
than 1.5 eV, which is roughly the difference between the bandwidth and the ionization po- 
tential of the clusters. Thus, the kinetic energies of the Auger electrons concentrate in a 
narrow energy range. Comparison with the experimental results by Gantefor et al || leads 
us to argue that the photoelectron signal which cannot be explained as coming from direct 
photoemission or thermionic effects is due to intraband Auger processes. Regarding a com- 
parison with experimental results by Reiners and Haberland ||], one can see in Fig. 2 that 
the Auger life-time of a vacancy in the ls-shell, ta(Is), of Na^ is approximately 100 ps, 
whereas T^(lp) ~ 10 ps. Since the time scale for evaporation is probably larger than these 
values, our calculations suggest an explanation for the fact that light induced evaporation 
in Nagx takes place within a photon-energy range which is smaller than the bandwidth. 

It is important to point out that our results are not sensitive to the particular form of the 
functional used for the exchange and correlation term. We have found neither qualitative 
nor appreciable quantitative differences by using the LDA functionals terms proposed by 
Wigner ||10|| , Gunnarson-Lundqvist |TI| , and the parametrizations by Perdew-Zunger fl2] 



and Vosko et al [jnj to the Monte-Carlo calculations of Ceperley and Alder JT4 
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We have also performed a Hartree-Fock (HF) calculation of the energy levels In 
contrast to the results obtained using LDA, the HF-treatment of the exchange term yields 
that Auger-transitions for neutral and positively charged clusters are energetically allowed. 
However, the number of possible transitions in Najv is much smaller than in Na^, and for 
Na^f there are just a few Auger-processes. Thus, the HF-calculation confirms the general 
trends obtained using LDA. 



IV. SUMMARY AND OUTLOOK 

We have calculated the intraband Auger-decay probabilities of spherical sodium clusters. 
Our results suggest that, for certain cluster sizes, there could be a competition between non- 
radiative electronic and atomic deexcitation channels. We found that the Auger-probabilities 
of small metal clusters are remarkably different from the atomic ones. In view of the re- 
sults presented in this paper it remains an interesting problem to study intraband Auger- 
probabilities in nanostructures and films and to compare them with the case of small clusters. 
Also the study of spin- dependent Auger-processes in small magnetic clusters appears as a 
possible interesting extension of this work. 
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FIGURES 

FIG. 1. Schematic illustration of the Auger emission in a Najy cluster with a vacancy in a deep 
electronic shell. Note that the Kohn-Sham effective potential V e ff(r) shows a barrier, in contrast to 
the case of neutral or positively charged clusters. The Auger probability for this particular process 
is given by the exchange matrix element (k, n"£"\V\n£, n'i') (see text). The direct Auger transition 
is obtained by exchanging the indices n'i' and nl. 

FIG. 2. Size dependence of the Auger life-time (in picoseconds) for a hole in the a) first shell 
(Is), b) second shell (lp) of spherical Na^ clusters. 

FIG. 3. Auger life-time for holes on different shells of Naj~ 37 . The dashed line refers to 
calculations taking into account relaxation effects for the final state. 

FIG. 4. Auger electron distribution P(s) as a function of the electron kinetic energy for a) 
NaJ" 37 , and b) Na^ 53 . 
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